Antiferromagnetic superexchange interactions in LaOFeAs 
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From first-principles calculations, we have studied the electronic and magnetic structures of the 
ground state of LaOFeAs. The Fe spins are found to be collinear antiferromagnetic ordered, re- 
sulting from the interplay between the strong nearest and next-nearest neighbor superexchange 
antiferromagnetic interactions. The structure transition observed by neutron scattering is shown 
to be magnetically driven. Our study suggests that the antiferromagnetic fluctuation plays an im- 
portant role in the Fe-based superconductors. This sheds light on the understanding of the pairing 
mechanism in these materials. 
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Recently an iron-based material LaOFeAs was re- 
ported to show superconductivity with a transition tem- 
perature T c ~ 26K by partial substitution of O with F 
atoms [1]. Soon after, other families of Fe-As oxyarsenides 
with La replaced by Sm[2[, Ce[3[, Pr[4[ and other rare 
earth elements were found superconducting with T c more 
than 50K. Like cuprates, these iron arsenides have a lay- 
ered structure. The superconducting pairing is believed 
to happen in the iron-based FeAs layers. The high tran- 
sition temperature and the preliminary band structure 
calculation suggests that the superconductivity in these 
Fe-arsenide superconductors is not mediated by electron- 
phonon interaction. It is commonly believed that the 
understanding of electronic structures of the parent com- 
pound LaOFeAs is the key to determine the underlying 
mechanism to make it superconducting upon doping. 

The early band structure calculations suggested that 
the pure LaOFeAs compound is a nonmagnetic metal but 
with strong ferromagnetic or antiferromangtic (AFM) in- 
stability p,@, 0]. Later, it was found that the antiferro- 
magnetically ordered state [1,0] has a lower energy than 
the nonmagnetic one, probably due to the Fermi surface 
nesting^. Dong et al.[10] predicted that the AFM state 
should form a collinear-striped structure by breaking the 
rotational symmetry. This collinear ordered AFM state 
has indeed been observed by the neutron scattering ex- 
periment 11, [l2[. Furthermore, the neutron scattering 
measurement found that there is a structure transition 
with a monoclinic lattice distortion at ~ 150K and the 
collinear order is formed about 15~20K below this tran- 
sition temperature. Without this distortion, the square 
AFM order induced purely by the Fermi surface nesting 
is expected to be more stable since there are two orthogo- 
nal but equivalent nesting directions (7r, tt) and (7r, — 7r), 
which can lower the energy of the ground state by keeping 
its rotational symmetry [8]. 

In this paper, we report the theoretical result on the 
electronic and magnetic structures of the ground state 
of LaOFeAs obtained from first-principles band struc- 



ture calculations. We find that there are strong near- 
est and next-nearest neighbor superexchange interac- 
tions in this material (similar conclusion was obtained 
by Yildirim[13[). The nearest and next nearest neighbor 
superexchange interactions have almost the same ampli- 
tude within error of calculation. Their competition af- 
fects strongly the electronic structure of the ground state. 
This drives a small monoclinic lattice distortion and a 
collinear ordering of Fe spins, as observed by neutron 
scattering. 

In our calculations the plane wave basis method was 
used{l4l]. We adopted the local (spin) density approx- 
imation and the generalized gradient approximation of 
Perdew-Burke-Ernzerhof[15j for the exchange-correlation 
potentials. The ultrasoft pseudopotentials [16[ were used 
to model the electron- ion interactions. After the full con- 
vergence test, the kinetic energy cut-off and the charge 
density cut-off of the plane wave basis were chosen to be 
600eV and 4800eV, respectively. The Gaussian broad- 
ening technique was used and a mesh of 16 x 16 x 8 
k-points were sampled for the Brillouin-zone integra- 
tion. LaOFeAs has a tetragonal layered structure with 
PA/nmm symmetry. A crystal unit cell consists of eight 
atoms with alternating FeAs and LaO layers along the c 
axis. In the calculation, the internal atomic coordinates 
within the cell were determined by the energy minimiza- 
tion. 

The previous band structure calculations @, @, 0, H, 
used an a x a x c crystal unit cell as the working cell, in 
which two Fe, As, La, and O atoms were included. To 
explore the magnetic structure, in particular the collinear 
AFM state of LaOFeAs, here we use a \[2a x \J~2a x c 
unit cell (Fig. 1). In order to determine the values 
of the nearest and next-nearest neighbor coupling con- 
stants of spin-spin interaction, J\ and J 2 (see Fig. 1), 
we have evaluated the minimal energies of four different 
magnetic states of Fe ions with constraints imposed if 
not stable. These four states have nonmagnetic, ferro- 
magnetic, square AFM, and collinear AFM orders, re- 
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FIG. 2: (Color online) Charge density distribution of 
LaOFeAs in the (001) plane crossing Fe-Fe atoms (a) and 
in the (110) plane crossing Fe-As-Fe atoms (b). 



FIG. 1: (Color online) Schematic top view of the FeAs layer 
in LaOFeAs. The small dashed square is an a x a unit cell 
while the large dashed square is a \[2a x \[2a unit cell. The 
collinear ordered Fe spins in the ground state are shown by 
red arrows. 



spectively. If the energy of the nonmagnetic state of 
LaOFeAs is set to zero, we find that the energies of the 
ferromagnetic, square AFM, and colliear AFM states are 
(0.0905, -0.010875, -0.21475) eV/Fe, respectively. Thus 
the ground state is a collinear-ordered AFM state, in 
agreement with the experimental observation [ill. [l2| . 

The magnetic moment around each Fe atom is found to 
be about 2.2 ~ 2.6 /ig, varying weakly in the above three 
magnetically ordered states. This suggests that the spin 
of Fe ions is between 1 and 3/2. The magnetic moment 
obtained from the neutron scattering model is about 6 
times smaller than this result. This is probably because 
the correlated effect, especially the strong competition 
between different AFM states, has not been fully included 
in the density functional theory calculation. 

To quantify the AFM interactions in this material, we 
assume that these energy differences are predominantly 
contributed from the interactions between the Fe spins 
which can be modeled by the following frustrated Heisen- 
berg model with the nearest and next-nearest neighbor 
couplings Ji and J2 



h — j\ Si • Sj + J2 Si-Sj, 



(i) 



whereas (ij) and <C ij ^> denote the summation over 
the nearest and next-nearest neighbors, respectively. 
From the calculated energy data, we find that J\ ~ 
0M98eV/S 2 and J 2 - 0.05l0eV/S 2 (the detail how these 
parameters are determined is given in the Appendix). If 
the spin of each Fe ion 5=1, then J\ ~ 0.0498eF and 
J 2 - 0.0510eV. 

The above result indicates that there are competing 
AFM interactions between Fe spins. In particular, the 
AFM coupling between two next-nearest neighboring Fe 
spins is very strong. This will frustrate the spin Neel 
structure and give rise to a collinear ordered AFM ground 
state. 
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FIG. 3: (Color online) Total and orbital-resolved partial den- 
sity of states (spin-up part) of LaOFeAs in the stripe-ordered 
antiferromagnetic state. The Fermi energy is set to zero. 



To explore the origin of these AFM interactions, we 
have calculated the charge distribution around Fe and As 
ions. The result (see Fig. [2^ and [2)3) shows that there 
is almost no charge distribution between two diagonal 
Fe atoms, but there is a strong bonding between Fe and 
As ions. This indicates that the AFM coupling J2 is 
induced by the superexchange bridged by As ions. This 
superexchange is AFM because the intermediated state 
associated with the virtual hopping bridged by As ions 
is a spin singlet. 

The charge distribution between two nearest Fe ions 
is finite. Thus there is a direct exchange interaction be- 
tween two neighboring Fe spins. Since there is a strong 
Hund's coupling between the spins of 3d electrons within 
each Fe ion, the direct exchange interaction is found to 
be ferromagnetic when the distance of two Fe atoms is 
between 2.4Aand 2. 85 A. However, the overall magnetic 
coupling J\ between the two nearest Fe spins in LaOFeAs 
is antiferromagnetic. Thus J\ is also dominated by the 
superexchange interaction bridged by As 4p orbit als. 

In the collinear AFM phase, we find that a small struc- 
tural relaxation, for which the lattice is slightly expanded 
along the antiferromagnic ordering direction (y-axis in 
Fig. 1) and shrunk along the ferromagnetic ordering di- 
rection (x-axis in Fig. 1), can further lower the ground 
state energy. This changes the angle between two prin- 
cipal axes in ab-plane, 7, from 90° to 90.47°. The corre- 
sponding energy gain is 7meV. However, we find that this 
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FIG. 4: (Color online) Total and projected density of states 
at the five Fe-3<i orbit als around one Fe atom. The Fermi 
energy is set to zero. 



small lattice distortion affects weakly the band structure 
and the Fe moments. 

Fig. [3] shows the total and projected density of states 
of LaOFeAs at the collinear AFM phase. In comparison 
with the nonmagnetic phase, we find that most of the 
states around the Fermi level are gapped by the collinear 
AFM order. This suppresses the total carrier density by 
more than two orders of magnitude. The strong suppres- 
sion is consistent with the Hall coefficient measurement 
which shows that the absolute value of the Hall coefficient 
is enhanced by more than 150 times in the AFM phase 
at 4K versus the nonmagnetic phase above 150K. Fur- 
thermore, we find that the density of states at the Fermi 
level is also suppressed compared with the nonmagnetic 
state[8|- However, it is not suppressed as strongly as for 
the total carrier density. The corresponding electronic 
specific heats are evaluated as 0.65mJ/(K 2 *mol) (stripe- 
ordered AFM) and 4.28mJ/(K 2 *mol) (nonmagnetic), re- 
spectively. This is also consistent with the specific heat 
measurement. 

We have also calculated the band structure of LaOFeAs 
with 5% F- or 5% Sr-doping by taking the virtual crystal 
approximation in the collinear AFM phase. We find that 
the overall band structure is hardly changed by 5% elec- 
tron or hole doping. The Fe moment is also unchanged. 
Only the Fermi energy moves up or down with electron 
or hole doping. However, as shown in Fig. [5] the Fermi 
surface changes dramatically. 

By projecting the density of states onto the five 3d 
orbitals of Fe (Fig. [4j), we find that the five up-spin or- 
bitals are almost completely filled and the five down-spin 
orbitals are only partially filled. However, the down-spin 
electrons are nearly uniformly distributed in these five 
3d orbitals. This indicates that the crystal field split- 
ting imposed by As atoms is relatively small and the Fe 
3d-orbitals hybridize strongly with each other. As the 
Hund rule coupling is strong, this would lead to a large 




FIG. 5: (Color online) (a) The electronic band structure 
of LaOFeAs in the collinear-ordered antiferromagnetic state 
with the in-plane angle 7 = 90.47° (see Fig. 1 for the def- 
inition of 7). The Fermi energy is set to zero. E' F and Ep 
correspond to 5% F- and 5% Sr-doping cases, respectively, 
(b) The Brillouin zone, (c) The Fermi surface of the un- 
doped compound: a hole- type cylinder along TZ and two 
electron- type pockets between F and X. FX corresponds to 
the parallel- aligned moment line and FX corresponds to the 
antiparallel-aligned moment line, (d) The Fermi surface of 
the 5% F-doped compound, (e) The Fermi surface of the 5% 
Sr-doped compound. 



magnetic moment formed around each Fe atom, as found 
in our calculations. The frustration between the J\ and 
J2 terms will suppress strongly the AFM ordering at the 
two Fe-sublattices, each connected only by the J 2 terms. 
This, together with the quantum fluctuation, will reduce 
strongly the average magnetic moment around each Fe 
measured by experiments. 

Fig. [5] shows the band structure and the Fermi sur- 
faces of electrons in the collinear AFM state. Unlike in 
the nonmagnetic state, there are now only three Fermi 
surface sheets in undoped case, one small hole cylinder 
along TZ and two small electron pockets formed between 
T and X, crossing the Fermi level. From the volumes 
enclosed by these Fermi surfaces, we find that the hole 
carrier density is about 1.64 x 10 19 /cm 3 and the electron 
carrier density is about 0.94 x 10 19 /cm 3 . Both decrease 
by more than two orders of magnitude in comparison 
with the nonmagnetic or square- AFM states [8[. Upon 
F(Sr)-doping, the electron (hole) Fermi surface sheets 
expand while the hole (electron) Fermi surface sheets 
shrink. With 5% F(Sr)-doping, the whole Fermi sur- 
face becomes electron(hole)-typed and the correspond- 
ing electron(hole) carrier density is 6.31 x 10 20 /cm 3 
(7.01 x 10 20 /cm 3 ), increasing by about 25 times com- 
pared with total carrier density in the undoped case. 

The above discussion shows that there are strong near- 
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est and next nearest neighbor superexchange interactions 
in LaOFeAs. The interplay between these AFM interac- 
tions can affect strongly the magnetic structure of the 
ground state. Upon doping, the AFM ordering will be 
suppressed. However, we believe that the remanent AFM 
fluctuation will survive, similar as in cuprate supercon- 
ductors. Thus the effective low energy model for describ- 
ing these Fe-based superconductors, no matter whether 
it is a single- or multi-band Hamiltonian, should include 
the frustrated Heisenberg terms defined by Eq. ([Tj). 

In the Fe-based superconductors, the magnetic fluc- 
tuation, induced by either the AFM superexchange in- 
teractions or the on-site Hund's rule coupling, can be 
responsible for the superconducting pairing. The for- 
mer interaction favors a spin singlet pairing, while the 
latter favors a spin triplet pairing. However, the super- 
conductivity induced by the Hund's rule coupling would 
generally involve the interband pairing, which is limited 
by the available phase space if the total momentum of 
Cooper pair is zero. This would suggest that the spin 
triplet pairing is not energetically favorable in a system 
with strong AFM fluctuations. Moreover, from the study 
of high-T c cuprate superconductivity, we know that the 
strong next-nearest neighbor AFM interaction favors an 
extended s or d xy -wdive pairing. Therefore, we believe 
that the leading pairing instability will be in spin singlet 
channel, if the superconductivity is driven by the AFM 
fluctuation. However, the competing nearest neighbor 
AFM interaction may introduce a small new component 
with different symmetry, for example a d x i _ y 2 -wave gap, 
to the pairing function. Thus the resulting gap parame- 
ter will generally be a superposition of two components 
with different symmetries, for example an extended s plus 



pairing state. 



In conclusion, we have presented first-principles calcu- 
lations of the electronic structure of LaOFeAs. We find 
that there are strong antiferromagnetic nearest and next- 
nearest neighbor superexchange interactions, bridged by 
As 4p orbit als. The next nearest neighbor antiferromag- 
netic coupling is comparable to the nearest neighbor one. 
This gives rise to the collinear AFM ordering of Fe spins 
in the ground state as observed by neutron scattering. 
The existence of strong antiferromagnetic fluctuations in 
the Fe-based superconductors bears a strong analogy to 
the high-T c cuprates. This suggests that the supercon- 
ductivity in these two different kinds of high-T c materials 
may have a common origin. 

We wish to thank N.L. Wang, J.L. Luo, P. Dai, Y.H. Su 
and H.G. Luo for fruitful discussions. This work is par- 
tially supported by National Natural Science Foundation 
of China and by National Program for Basic Research 
of MOST, China. While this paper was being final- 
ized, we learnt of the work of S. Ishibashi and coworkers 
(arXiv:0804.2963). They obtained a similar band struc- 
ture. 

Note added in revision: In a 2 x 2 Fe-Fe lattice of the 
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FIG. 6: The three magnetic configurations and the magnetic 
bond energies denned in the text. 



J1-J2 model with periodic boundary conditions, the J\ 
and J 2 terms are overcounted by a factor of 2 and 4, 
respectively. This overcounting has not been considered 
in our previous estimation of J\ and J2, which leads to 
an overestimation of the values of these parameters as 
well as the ratio J2 / J\ . This error has been corrected in 
this revised version. 



APPENDIX: DETERMINATION OF J x AND J 2 

To determine the value of J\ , one need to first evaluate 
the energy of a pair of nearest Fe-Fe moments in parallel 
(Ef,i) and anti-parallel (Ea,i) alignments with respect to 
a non-magnetic reference state, respectively. Then from 
their difference, one can determine the value of J\ by the 
following formula 



Ji = (E F>1 - E AA )/(2S 2 ), 



(2) 



where S is the spin of each Fe ion. It should be empha- 
sized that Ep,i is not necessary to be equal to — Ea,i 
since the energy of the reference state may not be lo- 
cated exactly at the middle of the energy between the 
ferromagnetic and antiferromagnetic states. This energy 
lineup between the nonmagnetic and any other magnetic 
state needs self-consistent total energy calculations to de- 
termine, that is what we have done. Thus Ep,i and Ea,i 
should be determined independent. Similarly, J 2 can be 
determined from the difference between the energy of a 
pair of next-nearest Fe-Fe moments in the parallel (Ep,2) 
and anti-parallel (Ea,2) alignments: 



J2 = (E F , 2 - E Aa )/(2S 2 ). 



(3) 



To determine the values of Ep,i — Ea,i, Ep,2 — Ea,2, 
we have calculated the total energies of the ferromagnetic 
(Efm), square antiferromagnetic (Eaf), an d collinear 
antiferromagnetic (Ecol) states, respectively. The spin 
configurations of these three states are shown in Fig. 
[6j The corresponding energy differences with respect to 
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indeed due to the competition of superexchange interac- 
tions. 



The above estimation indicates that J\ ~ J2 within 
the error of calculation for LaOFeAs. In this param- 
eter range, as shown by Yao and Carlson [17] with the 
spin wave approximation, the J\ superexchange interac- 
tion term will compete strongly with the J2-term. This 
results in a strong reduction of the magnetic moment of 
Fe, which would naturally explain why the observed mag- 
netic moment (~ 0.36 \±b) is significantly smaller than 
that obtained from the density functional calculations 
for LaOFeAs. 



FIG. 7: Variations of J\ and J2 with Efm — Enm- The 
dashed line denotes our calculated value. 

the nonmagnetic state (Enm) are (0.0905, -0.10875, - 
0.21475) eV/Fe. If these energy differences result mainly 
from the exchange interactions between the nearest or 
next-nearest Fe moments, then we obtain the following 
equations 

E F m-E N m = 2E Fjl +2E Fj2 = 0-0905eV, 
E AF - E NM = 2E AA + 2E F , 2 = -0.10875eF, 
Ecoi - E NM = E Ffl + E A ,i + 2E A , 2 = -0.21475eV. 

From them, we further find that 

E FA - E A ,i = 0.0996el/, 
E Fy2 -E Ay2 = 0.1028eV. 

Thus the values of J\ and J2 are 

Ji = 0.0498eF/S 2 , (4) 
J 2 = 0.050W/S 2 . (5) 

The energy of the current ferromagnetic state is less ac- 
curately determined since this state is not a stable state. 
The error in E F m — Enm will give rise to the error in 
J 1 and J 2- Fig. [3 shows how J\ and J 2 change with 
E F m — Enm- As we see, we find that J2 > Ji/2 even 
we assume the deviation of E F m — Enm from our calcu- 
lated value is as big as 0.1 eV. Thus we believe that the 
collinear antiferromagnetic order observed in LaOFeAs is 
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